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Cross  section  area  of  exit  of  nozzle,  in^ 

2 

Cross  section  area  of  ejector,  in 
Diameter  of  exit  of  nozzle,  in 
Diameter  of  ejector,  in 
Slope  of  ^  vs  %  curve 
Length  of  gap,  in 
Length  of  ejector,  in 

Inverse  of  exponent  related  to  velocity  profile  at  outlet  of 
ejector,  defined  in  Equation  2 

Exponent  related  to  velocity  profile,  defined  in  Equation  36 

2 

Local  pressure,  lb/ in 

2 

Pressure  of  the  atmosphere,  lb/in 

Pressure  at  exit  of  nozzle,  lb/in^ 

Pressure  at  initial  plane,  lb/in^ 

Pressure  at  outlet  of  ejector,  lb/in^ 

Reynolds  Number, 

Radius  of  exit  of  nozzle,  in 
Radius  of  ejector,  in 

Distance  measured  in  radial  direction  from  ejector  centerline, 
in 

Radius  of  the  laminar  core,  in 

Thrust  force  of  jet  ejector  with  no  gap,  lb 

Thrust  force  of  jet  ejector  with  gap,  lb 


Thrust  force  on  jet  flow,  lb 

Axial  velocity,  ft/sec 

Centerline  axial  velocity,  ft/sec 

Axial  velocity  at  the  exit  of  nozzle,  ft/sec 

Secondar,  or  induced  flow  velocity  at  initial  plane,  ft/sec 

Velocity  at  the  outlet  of  ejector  and  its  centerline  value,  ft/sec 

Potential  velocity  beyond  the  edge  of  the  jet,  ft/sec 

Radial  velocity,  ft/sec 

Distance  in  radial  direction  from  edge  of  the  laminar  core,  in 

Distance  in  radial  direction  from  the  edge  of  the  laminar  core 
to  the  edge  of  jet,  in 

Axial  distance  from  the  initial  plane,  in 


Length  of  the  laminar  core,  in 

Length  of  the  modified  free  mixing  flow  regime,  in 

Distance  in  the  radial  direction  from  the  centerline  of  the 
ejector  to  the  edge  of  the  jet,  in 

Radius  of  the  jet  at  x  =  x^,  in 
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Shear  stress,  lb/in^ 

Thickness  of  boundary  layer  of  the  walls  of  ejector,  in 
Kinetic  viscosity  of  air,  ft'vsec 


INTRODUCTION 


A  jet  ejectro,  or  jet  pump  as  it  is  sometimes  called, is  a  device  in 
which  a  fluid  is  emitted  from  an  orifice  or  nozzle  exit  into  a  hollow, 
generally  cylindrical,  body.  Due  to  the  shear  stresses  between  the 
moving  fluid  and  the  ambient  still  fluid,  part  of  the  kinetic  energy  of 
the  moving  is  transferred  to  the  still  fluid.  In  the  case  of  a  constant- 
pressure  mixing  tube,  which  can  be  made  through  careful  design, 

(Reference  1),  this  mixing  process  is  the  primary  cause  of  secondary  or 
induced  flow. 

However,  in  the  case  of  a  constant-area  mixing  tube,  a  low-pressure 
field  is  created  by  the  jet  flow.  This  difference  in  pressure  between 
the  exit  plane  of  the  nozzle  and  the  ambient  fluid  upstream  causes  an 
entrainment  of  the  ambient  fluid  or  a  secondary  flow.  Due  to  the  trans¬ 
fer  of  energy,  the  secondary  flow  is  either  accelerated  or  pumped  to  a 
higher  pressure.  The  jet  flow  and  its  induced  flow  ate  not  necessarily 
of  the  same  fluid;  either  of  them  can  be  air,  water,  steam,  or  the 
exhausted  gas  of  a  jet  engine.  The  ejector  is  not  necessarily  circular, 
but  can  be  of  any  shape,  as  required  by  the  physical  installation. 

The  principle  of  a  confined  jet  flow  finds  many  different  applica¬ 
tions  in  industry,  especially  the  aerospace  industry.  An  ejector  can 
be  used  as  a  pumo,  blower-augmenter ,  and  noise-reducer.  (References  2, 

3,  4,  5,  6,  7,  8,  9,  10,  11,  12).  In  case  the  jet  ejector  is  operated 
with  no  change  in  pressure  it  does  not  provide  a  pumping  effect.  How¬ 
ever,  it  does  provide  a  thrust  increase  for  the  installation. 

A  typical  constant-area  jet  ejector  arrangement  is  shown  schemat¬ 
ically  in  Figure  1.  With  this  arrangement,  the  primary  jet  flow  is 
emitted  from  the  exit  of  the  nozzle  into  the  ejector,  and  the  velocity 
is  almost  uniformly  distributed  throughout  the  cross  section  area  of  the 
exit  of  the  nozzle.  A  subatmospheric  pressure  field  is  created,  which 
consequently  causes  a  secondary  flow  in  the  rest  of  the  cross  section 
area  of  the  ejector  at  the  initial  plane.  If  the  boundary  layer  growth 
along  the  contour  of  the.  faired  inlet  of  the  ejector  is  neglected,  the 
secondary  velocity  distribution  can  be  considered  to  be  constant. 

The  velocity  of  the  secondary  flow  is  less  than  that  of  the  primary 
flow.  As  a  consequence ,  there  is  a  near  discontinuity  in  the  velocity 
profile  across  the  ejector  in  the  initial  plane.  Downstream  of  this 
plane,  a  process  of  turbulent  mixing  tends  to  smooth  out  this  near  dis¬ 
continuity  so  that,  if  the  ejector  is  long  enough,  the  mixed  flow  velocity 
profile  becomes  that  of  fully  developed  turbulent  pipe  flow. 

Since  the  static  pressure  increases  downstream  of  the  initial  plane 
until,  at  the  outlet  of  the  ejector,  it  becomes  equal  to  that  of  the 
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surrounding  atmosphere,  the  static  pressure  increase  per  unit  length  of 
the  ejector  depends  upon  the  total  length  of  the  ejector.  The  mass  rate 
flow  remains  constant  throughout  the  length  of  the  ejector,  but  total 
momentum  at  the  outlet  of  the  ejector  is  less  than  that  at  the  initial 
plane  due  to  losses  which  result  fi^m  skin  friction. 

The  maximum  augmentation  factor  for  a  straight  ejector  is  2.0  as 
calculated  theoretically  by  von  Karman.  (Reference  6).  He  assumed  that 
the  velocity  distributed  was  uniform  at  the  outlet  of  the  ejector,  and 
that  skin  friction  losses  were  negligible.  A  similar  calculation  has 
been  performed  in  this  report  using  somewhat  different  assumptions. 

A  knowledge  of  the  free-mixing  process  and  pipe  flows  aids  in 
understanding  the  confined-mixing  process.  Information  concerning  those 
phenomena  can  be  found  in  many  works.  (References  13,  14,  15,  16,  17, 

18,  19,  20,  21).  The  equations  of  motion  for  the  flow  field  of  a 
constant-pressure  mixing  tube  is  similar  to  that  of  a  free  jet  flow, 
however,  the  boundary  conditions  are  different.  The  axial  pressure 
gradient  term  in  the  equation  of  motion  for  a  constant-area  mixing  tube 
and  the  coexistence  of  laminar  and  turbulent  flows  make  the  flow  field 
difficult  for  mathematical  analysis. 

For  comparison  purposes  a  satisfactory  solution  of  the  problem  of 
a  constant-area  mixing  process  can  be  obtained  by  approximate  solution 
of  the  Naview-Stokes  equations.  However,  this  approach  requires  such  a 
great  amount  of  work  that  it  becomes  questionable  whether  the  approxi¬ 
mate  solution  justifies  the  effort.  As  a  consequence,  this  report  adopts 
an  approach  in  which  the  author  uses  existing  theories  and  empirical 
formulae  to  solve  the  problem. 
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EXPERIMENTAL  APPARATUS  AND  TECHNIQUES 


For  the  sake  of  convenience,  tTo  sets  of  experimental  apparatus 
were  used  -  one  for  measuring  the  thrust  forces  and  one  for  measuring 
the  velocity  distributions.  Two  convergent  nozzles  and  three  constant- 
area  ejectors  were  used  in  both  measurements;  with  one  constant-area 
ejector  of  relatively  long  length  being  used  in  velocity  distribution 
measurements  only.  Compressed  air  used  in  both  measurements  was 
supplied  by  the  equipment  shown  schematically  in  Figure  2.  Thrust 
forces  were  measured  with  the  device  shown  in  Figure  3a. 

In  order  to  change  the  relative  location  of  the  exit  of  the  nozzle, 
brass  tubes  were  used  as  extensions  of  the  nozzle  exit.  (Figure  4). 
Their  dimensions  are  shown  in  Table  1. 

Ejectors  of  various  lengths  were  attached  to  the  primary  nozzle  as 
shown  in  Figure  3b.  Care  was  taken  to  insure  that  the  ejector  and 
primary  nozzle  were  aligned  axially.  The  dimensions  of  the  ejectors 
are  shown  in  Table  2. 

By  using  these  nozzles  and  ejectors,  the  area  ratio  of  the  jet 
ejector  could  be  varied  from  4.0  to  105.5. 

The  primary  nozzle  flow  and  the  secondary  entrained  flow  could  be 
calculated  from  velocity  profiles  measured  at  various  locations  in  the 
ejector  by  means  of  a  pitot-static  tube.  This  tube  was  attached  to  a 
micrometer-type  mechanism  which  allowed  it  to  be  positioned  at  any 
desired  location  along  the  diameter  of  the  ejector.  (Figure  5).  Velo¬ 
cities  were  measured  by  means  of  a  U-tube  water  manometer. 

Profiles  were  measured  at  the  distances  downstream  of  the  initial 
plane  as  shown  in  Table  3. 

All  velocity  profiles  and  static  pressure  were  measured  when  the 
nozzle  and  ejector  were  positioned  horizontally.  (Figure  6).  All 
thrust  force  measurements  were  made  with  the  apparatus  positioned 
vertically. 

Thrust  forces  of  the  nozzles  were  firsc  measured  at  various  gauge 
pressures  as  shown  in  Figure  7.  These  measurements  :ere  hen  repeated 
with  the  ejector  attached  and  were  plotted  against  gauge  pressure. 
Typical  curves  are  shown  in  Figure  8.  At  an  arbitrary  gauge  pressure  of 
35  p.s.i.  the  augmentation  factor  was  obtained  by  dividing  the  thrust 
of  the  jet  with  ejector  by  the  thrust  of  the  jet  without  ejector.  The 
effect  of  the  length  ratio  on  augmentation  factor  is  shown  in  Figure  9. 


The  thrust  forces  of  the  jet  ejectors  at  different  gaps  were 
measured.  These  forces  were  divided  by  the  thrust  cf  the  jet  ejectors 
with  no  gap  in  order  to  obtain  the  gap  augmentation  factor  (pQ  ,  and 
were  similarily  divided  by  the  thrust  of  the  jet  in  order  to  obtain  the 
overall  augmentation  factor,  Qr-tto  and  QV  are  plotted  against 
gap  ratio  in  Figure  10.  Gap  is  the  distance  between  the  nozzle  exit 
plane  and  the  initial  plane  of  the  ejector.  The  initial  plane  is 
defined  as  the  plane  of  the  ejector  where  the  faired  inlet  of  the 
ejector  ends  and  the  constant-area  portion  of  the  ejector  begins.  Gap 
is  defined  as  negative  when  the  exit  of  the  nozzle  is  inserted  into  the 
ejector,  as  positive  when  the  exit  is  away  from  the  ejector,  and  as  zero 
when  the  nozzle  exit  plane  coincides  with  the  initial  plane  of  the 
ejector. 


THEORETICAL  CONSIDERATIONS 


Several  assumptions  are  made  in  the  theoretical  considerations. 

They  are : 

1.  The  flow  is  incompressible,  hence  the  static  pressure  is 
constant  throughout  the  initial  plane  within  the  ejector 
and  the  pressure  at  the  outlet  of  the  ejector  is  equal  to 
that  of  the  surrounding  atmosphere. 

2.  The  velocity  at  the  exit  of  the  nozzle  is  uniformly  distributed, 

i.e.,  the  velocity  profile  is  constant 

3.  The  secondary  flow  velocity  is  uniformly  distributed  at  the 
initial  plane 

4.  There  is  no  pressure  change  in  the  radial  direction. 

5.  The  flow  field  is  of  an  axisymmetric  case. 


A,  Thrust  Augmentation  Factor  and  Velocity  Ratios 

A  method  similar  to  that  used  by  von  Karman  (Reference  6)  was  used 
in  calculating  the  thrust  augmentation  factor  and  the  velocity  ratios. 
It  was  also  assumed  that 


1.  The  velocity  distribution  at  the  exit  of  the  ejector  follows 
the  relation;  ..  I 

u*  li— 
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2.  Skin  friction  losses  are  negligible, 

3.  The  walls  of  the  nozzLe  are  infinitely  thin  at  the  initial  plane 

4.  The  velocities  at  the  nozzle  exit  remain  unchanged  by  the  instal 
lation  of  the  ejector. 


The  equation  of  continuity  may  be  written  as r 

PAUc+f(Ae-A)  U,-  P  f^irUaZ-Jz 


After  dividing  throughout  by  one  may  write  the  equation  as 

/  U.a.\  f  Z 

'o 


i+ <,*-/;>•< = ,2^  r'  d  <£•> 
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Using  the  assumption  that 


CtWfe) 


Z  \w 


(2) 


equation  1  may  be  rewritten  as 


2f'rr\ 


I  +  ^-1)^  ,  )X.<z  . 

The  momentum  equation  may  be  written  as 

AfPUo**  ?o)  +  (Ae~A)((>iJl:L+1Pl)=J*e((3\lz+i>2)ZTTZ:clZ; 

and  according  to  Bernoulli's  theory 


(3) 


where 


frfo,  Pz=r* 
pa.*?’  +  *PU* 


The  momentum  equation  may  be  rewritten  as 

I  +  i°<r  (&e)ol  (We). 

Integrating  the  term  on  the  right-hand  side  results  in 

i  +  z*?  (a~2) ~  i 

Solving  the  velocity  ratios  from  equations  3  and  5  for  m  =4, 

„  _  -Okffl+M  ¥7^  W^^kE^i 

”  a  —  ^  -  - -  —  . . 


(4) 


(3) 


and 
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li  vl.  is  much  greater  than  1.0,  equations  6  and  7  may  be  simplified  and 
rewritten  as 


(9) 


The  thrust  augmentation  factor,  u  ,  is  defined  by  the  formula 

in  —  Jo 

Ul  ~  PAUo* 

Carrying  out  the  calculation,  one  has  for  the  thrust  augmentation  factor, 

(10) 

It  is  somewhat  arbitrary  to  define  the  thrust  augmentation  factor 
as  the  ratio  of  the  actual  thrust  to  the  product  of  the  mass  and  the 
velocity  of  the  jet  instead  of  to  its  total  impulse-  However,  in 
general,  the  difference  is  not  large. 

In  Figure  12  the  thrust  augmentation  factor,  (Jf  ,  is  shown. 

Velocity  ratios  °<  i  and  are  plotted  in  Figure  13,  showing  compari¬ 

son  witt.  the  experimental  measurements. 

The  pressure  reduction  may  be  obtained  by  using  equation  4 


This  equation  has  been  plotted  in  Figure  14. 


B.  Structure  sL£h&  Flow  Field  in  the  Inlet  Length  of  the  Elector 

The  structure  of  the  flow  field  at  the  inlet  of  the  ejector  is  very 
complicated  (Figure  1).  As  the  compressed  air  is  emitted  from  the  nozzle 
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into  the  constant-area  ejector,  there  is  a  field  of  subatmospheric 
pressure  created.  Consequently,  a  secondary  flow  of  constant  velocity 
is  formed  at  the  inlet  of  the  ejector.  Disregarding  the  faired  portion  of 
the  inlet,  the  boundary  layers  at  the  wall  are  infinitely  thin  at  the 
initial  plane.  Normal  boundary  layer  growth  occurs  downstream  from  the 
initial  plane.  Disregarding  the  flow  of  the  jet,  the  boundary  layers 
keep  growing,  until  they  meet  at  a  point  on  the  axis  of  the  pipe.  Until 
this  happens,  there  is  a  core  of  fluid  practically  uninfluenced  by 
viscous  effects  in  which  the  total  head  may  be  considered  constant. 

This  problem  of  inlet  flow  of  a  circular  pipe  had  been  investigated  by 
Goldstein  (Reference  22),  Langhaar  (Reference  23),  Schiller  (Reference 
24),  and  Talbot  (Reference  25). 

Due  to  the  existence  of  the  jet  flow,  the  edge  of  the  wall  boundary 
layers  will  not  meet  at  a  point  on  the  axis  of  the  pipe;  instead  they 
join  the  boundary  of  the  jet. 

The  flow  field  in  the  inlet  portion  of  the  ejector,  disregarding 
the  primary  jet  flow,  becomes  similar  to  the  case  studied  by  Schiller. 
Schiller  suggested  that  the  wall  boundary  layer  thickness,  £  ,  is 

related  to  the  velocities,  and  U^,  by  the  following  relationship, 

S=2Re-R  ^ 

where  Re  is  the  radius  of  the  circular  tube,  U  is  the  constant  velocity 
of  the  potential  core  at  the  initial  plane  andXU  is  the  velocity  in  the 
potential  core  which  is  a  function  of  x.  p 


The  velocities  have  been  related  to  the  axial  coordinate  as 


Re-  „  RTCn.), 


~  ~  J _  r~> 

4 


where 


o  =  ELeUji 
K  ^r 


^  _  Uf> 

"l-uf 


is  the  Reynolds  Number,  and 

-  I . 


The  function 


■fC»v  has  been  given  by  Schiller  as 

"f  Cv\)r  C^5  +  </4+Z  2  ' 

-  SiVi  - 


(12) 


(13) 


(14) 
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and  has  been  simplified  by  Goldstein  to  the  following  form  which  has 
been  plotted  in  Figure  11, 


8  [ff  -  Wj 4+Z>\rZ)rC  ‘ 

-  Lk/4-zn\i  3%fzT  -  I  zrx. Li2ck.-  rrn  ~ 

5  vTFyv'V  "  io  Sin  3  +  3  \oJz  Sin  3J. 


(15) 


For  the  sake  of  simplicity,  one  can  approximate  the  curve  of 
equation  15  uy 


-f(yj)  =  0.00133  +  Q'Ob(oZ5Y£  +  O.OOZObrtf . 


(16) 


Equation  16  is  plotted  in  Figure  11  for  comparison  with  equation  15. 


Since  the  thickness  of  the  wall  boundary  layer,  £  ,  is  an  explicit 

function  of  Up  (equation  11),  Up  is  an  implicit  function  of  X 
(equation  12).  Consequently,  £  is  a  function  of  Xalso.  Vro 
plotted  against  the  nondimensional  variable  lK4r>  in  Figure  15. 


is 


The  variation  of  the  radius  of  the  mixing  zone  can  be  estimated  by 
a  straight  line.  Starting  with  the  equation  given  by  Relmbold  (Refer¬ 
ence  26),  for  a  mixing  tube  of  constant  pressure, 


1.72.3  -I-  0.4fc7(l-°<i)  +  0.a5Cl-*f,y 


(17) 


where  denotes  the  radius  of  the  mixing  zone  at  the  end  of  the  laminar 
cjre  of  the  jet  flow.  Hence,  _ (g 


where  xfc, 
core,  is 


as  given  by  Helmbold,  expressing  the  length  of  the  laminar 


X*  ..  IU> 


(18) 


Thus 


X 

-5(1— <•.) 


[l.723  +  O.4fc7O-«(1)-fO.25(l-<!f,'>1-0+'1. 

(19) 


By  equating  the  radius  of  the  mixing  zone  (Re  -  $  ),  one  can  find 
the  location  where  the  wall  boundary  layer  joins  the  boundary  of  the 
mixing  zone . 
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The  primary  flow,  except  in  the  turbulent  boundary  layer  at  the 
wall  of  the  nozzle,  emerges  as  a  laminar  jet  with  uniform  velocity 
which,  beginning  at  the  rim  of  the  nozzle,  becomes  increasingly  turbulent 
downstream  as  the  mixing  of  primary  and  secondary  flow  spreads  inward. 

The  length  of  the  laminar  core  has  been  investigated  by  many  authors. 

For  a  free  jet,  according  to  Paris  (Reference  15),  it  is  4.3  nozzle 
diameters,  which  compares  favorably  with  Keuthe's  4.44  nozzle  diameters 
(Reference  14),  and  Davies'  4.35  nozzle  diameters  (Reference  '7);  but 
for  confined  jet  flow,  it  depends  on  area  ratio. 

For  the  confined  jet  flow  of  the  present  case,  the  dimensions  of 
the  core  could  be  computed  by  the  following  formulae  derived  from  the 
results  of  Szablesaki's  calculation  by  Helmbold  (Reference  26),  valid 
for  the  case  of  a  confined  jet  with  constant  pressure  across  the  ejector. 
The  length  of  the  laminar  core  in  terms  of  nozzle  radius  is 


The  radius  of  the  mixing  zone  at  the  distance  x^x^  is 

^-=1.723  +  0.467  (I—0  +  0.25  (l-*,)1' . 


(18) 


A  comparison  of  the  lengths  of  the  laminar  core  is  presented  in 
Figure  16. 


C.  Velocity  Distribution  Near  the  Inlet  of  the  Ejector 


The  velocity  distribution  of  a  confined  jet  flow  has  been  investi¬ 
gated  by  many  authors,  such  as  Cruse  and  Tontini  (Reference  1),  Helmbold 
(References  26,  27,  28)  and  Squire  and  Trouncer  (Reference  29). 


Since  the  flow  is  almost  parallel  and  the  pressure  is  constant  over 
any  cross  section  area  of  the  ejector,  the  simplified  form  of  the 
boundary  layer  equations  is  justified  here.  The  equation  of  motion  is 


UL 


<3  CL 

ax 


4-  tr  ^  =  —  X  4^-  +- 


pr 


a  r 


(21) 


Following  the  apparent  viscosity  theory  of  Prandtl  (Reference  21), 
the  shearing  stress,  ,  may  be  expressed  as  ^ 


ar  j 


(72) 
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where  £  ,  the  so  called  virtual  kinematic  viscosity,  could  he  written  as 

€  y  >  iiy ^23) 

H  is  a  constant,  and  Y  is  the  radius  of  the  boundary  of  the  jet  before 
the  boundary  meets  the  edges  of  the  boundary  layers  of  the  walls  of  the 
ejector.  Graax>  in  the  present  report,  is  the  centerline  velocity  U  and 
umin  is  the  potential  velocity  Up  outside  the  boundary  of  the  spreading 
of  the  jet. 

Thus,  the  shearing  stress,  following  the  terminology  used  in  this 
report,  could  be  written  as 

(24) 


Substituting  the  expression  for  the  shearing  stress  into  the 
original  equation  of  motion  (equation  21),  one  has 


a 


<^U- 

3'X 


^  *  t>  =  - +  *  y ' £r  [r (fr)J25) 


It  seems  reasonable  to  assume  the  radial  velocity  component,  v, 
disappears  at  the  outer  potential  flow  zone.  Since  the  axial  velocity 
component  is  a  function  of  x  only,  the  equation  of  motion,  within  the 
potential  flow  zone,  can  be  simplified  to 


(26) 


Consequently,  the  equation  of  motion  can  be  rewritten  as 

u  +u-f^=  u-p 


(27) 


As  suggested  by  Helmbold,  the  axial  velocity  component,  u,  may  be 
written  as 


u.  -  U  f  'x)  *  U  ( x'1  -f  (  x  j  r)  ^ 


(28) 


where  a  u  ( x; 


u  60- Uf  (*X 


1 


Substituting  equation  Z8  into  equation  //,  one  nas 

+  AUf  i )+  Lr^-Ufr* 

~  U'p  0-fo  -h  yc  y  *L)Z(^r~  rf)} 


where 


p/  /  _ 

trr  -  §fS  . 


pi  miL  Si  ~ 

IT*  >  r  sr,+x  *Xj> 


If  the  law  of  similarity  holds, 
a  or  JL—1  r  <  I _  __ 


fr'  =  yr«i^, 


The  boundary  conditions  are  as  follows: 

at  r»0  1  f  *  I  ,  hence  -fx  =  $ 


fr-0 
\y  ~  0 


at  r  =  y6<)  ; 


f-O  fr=0 


hence  _  '  _  0 
+  X  ~  ^ 


As  a  consequence, 
at  r  -  0  , 


CtL.  V 

(U-p-l-AU)(Up'+  aU')=U-pO*>  +'Xv(*'&UI%Z‘o 

XyM£To(¥+ fWWi+j&0s8r+^*}-. 
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At 


or 


/  ~N 

r  *  y  v  kj  j 

OpVp1  =  0r>Up  +•  ycyCfi^U2-  frtJ 

■f rtr  = 


(33) 


Helmbold  had,  for  the  sake  of  simplicity,  suggested  the  function 

f  I -(-$-)*  J4. 


This  function  satisfies  the  boundary  conditions  of  equation  33  and  of 
equations  30  and  31. 


Faris  (Reference  15),  in  studying  the  free  jet  mixing,  used  the 
velocity  distribution  suggested  by  Cornish  (Reference  13): 


(35) 


where  n  is  an  empirical  function  of  the  axial  coordinate,  x  varies 
between  1.0  and  2.0,  and  w^  denotes  the  w  where =^. 


Equation  35  should  be  changed  so  that  it  can  be  relevant 
compared  with  the  measurements  of  the  present  investigation. 
35  becomes 


U-—  Ufff_ 

V-U-P 


X[i-COSn{\-[l~  (o.S^JzOj-  }"J  , 


and 

Equation 


(36) 


U-U-p  _  ±, 

where  w^  denotes  the  w  where  y__T  ~  %, 
plotted  in  Figure  17  for  various  values  of 
experimental  measurements. 


.  Equation  36  is 
n  and  is  compared  with 


D.  Velocity  Distribution  in  the  Aft  Part  of  the  Ejector 

If  the  ejector  is  long  enough,  the  flow  eventually  becomes  fully 
developed  turbulent  pipe  flow.  That  is,  after  the  mixing  process  of 
the  primary  and  secondary  flows  in  the  inlet  length  of  the  mixing  tube 
and  the  transition  processes,  a  fluid  particle  will  enter  a  flow  regime 
where  the  velocity  profile  is  unchanged  with  respect  to  the  axial  dis¬ 
tance.  This  velocity  profile  can  be  expressed  as  a  function  of  the  radial 
coordinate  ,  Uz  v  ^  yfc 

VQT/  Ue/  : 

3  (38) 
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where  B~Re—Y  ,  and  m  is  a  function  of  the  area  ratio,  jjl.  ,  and  the 
length  ratio,  L/D,  or 

Experimental  measurements  and  Equation  38  are  compared  in  Figure 

18.  The  variation  of  m  with  respect  to  is  plotted  in  Figure 

19.  M- 
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DISCUSSION  AND  CONCLUSIONS 


A.  Discussion: 


The  experimental  measurements  for  optimum  thrust  augmentation  factor, 
velocity  ratios,  and  pressure  reduction  ratios  for  the  initial  plane  at 
various  area  ratios  are  compared  with  theory  in  Figures  12,  13,  and  14. 

The  discrepencies  between  experimental  measurement  and  theory  are 
smaller  for  the  low  area  ratios  than  for  the  high  area  ratios.  (Figure 
12).  This  phenomenon  can  be  explained  as  follows:  frictional  losses 
were  neglected  in  deriving  the  formula,  but  in  reality,  skin  friction 
losses  were  quite  influential.  For  low  area  ratios,  a  shorter  ejector 
was  adequate  to  produce  the  optimum  thrust  augmentation  factor;  but  for 
high  area  ratios,  a  longer  ejector  was  needed.  Consequently,  frictional 
losses  were  larger  for  the  high  area  ratios  than  for  the  low  area  ratios. 

The  uptimum  thrust  augmentation  factor  obtained  in  the  present  investiga¬ 
tion  was  1.50  when^Jb  was  105.5. 

Experimental  results  are  compared  with  theory  for  the  velocity  ratios, 

,  and  <*2  ,  in  Figure  13.  The  comparison  between  theory  and 
experiments  for  the  pressure  variation  with  respect  to  area  ratio  is 
presented  in  Figure  14.  The  augmentation  factor  assumed  its  maximum 
value  when  the  gap  ratio  was  0.8,  as  can  be  seen  in  Figure  10. 

It  seems  reasonable  to  divide  the  flow  field  in  the  ejector  into 
three  regimes.  (Figure  1).  The  first  regime,  wnere  the  mixing  process 
of  the  free  jet  is  almost  retained,  may  be  denoted  as  a  modified  free- 
mixing  regime.  Like  the  ca:  e  of  a  free  jet,  a  laminar  core  is  attached 
to  the  exit  of  the  nozzle.  Unlike  the  case  of  a  free  jet,  a  potential 
flow  parallel  to  the  axis  of  the  ejector  surrounds  the  laminar  core.  A 
near  discontinuity  in  velocity  exists  at  the  rim  of  the  exit  of  the 
nozzle.  This  near  discontinuity  smoothes  cut  downstream  due  to  the 
viscous  shearing  stresses.  The  mixing  zone  "survives"  the  laminar  core 
and  the  potential  zone  downstream,  and  eventually  becomes  fully  developed 
turbulent  pipe  flow.  (Figure  1) ,  The  wall  boundary  layer  starts  at  the 
rim  of  the  inlet  of  the  ejector  and  continues  to  grow  downstream.  The 
edge  of  the  wall  boundary  layer  meets  the  boundary  of  the  mixing  zone 
at  a  location  which  can  be  computed  theoretically  or  measured  experimentally. 

Theoretically , -%Ja-  is  11.0  for/t  =  16.0,  15.0  ioxJJi  =  29.6,  18.0 
for/L  =  46.7,  and  32.6  ioc JA-  -  105. 5.  Experimentally,  this  location  was 
found  to  be  close  to  the  theoretical  location  (Figure  15).  •~Tq~  is  not 
only  a  function  of  the  area  ratio,  but  is  also  affected  by  the  length 
ratio,  L/De.  Schiller,  in  his  theory  of  inlet  length  flow  of  circular 
pipes,  did  not  specify  the  length  of  the  pipe;  however,  he  implied  that 
the  pipe  length  was  so  great  that  the  conditions  at  the  inlet  portion 
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was  no  longer  a  function  of  the  length  of  the  pipe.  The  ejectors  of  the 
present  case  were  not  long  enough  to  fulfill  the  conditions  in  which 
Schiller's  theory  would  be  entirely  applicable.  The  length  of  the 
ejector  from  the  initial  plane  up  to  the  location  where  the  edge  of  the 
wall  boundary  layer  joins  the  edge  of  the  jet  flow  spread  is  denoted  as 
the  modified  free-mixing  flow  regime. 

As  mentioned  above,  the  flow  field  eventually  becomes  fully  developed 
turbulent  pipe  flow  if  the  ejector  is  long  enough.  A  fully  developed 
turbulent  pipe  flow,  as  denoted  in  this  report,  has  a  velocity  profile 
which  may  be  expressed  by  a  fraction  power  formula, 

XT*  ~  V  ReJ  } 

where  m  is  a  positive  real  number,  z  is  the  distance  measured  from  the 
walls,  and  Re  is  the  radius  of  the  pipe. 

The  flow  field  in  the  ejector,  after  the  potential  flow  vanishes 
and  before  the  flow  field  becomes  fully  developed  turbulent  pipe  flow, 
is  denoted  as  the  transition  regime.  The  portion  of  the  ejector,  from 
the  location  where  the  velocity  profile  becomes  that  of  a  turbulent  pipe 
flow  to  the  end  of  the  ejector,  is  denoted  as  the  pipe  flow  regime. 
However,  since  the  velocity  profiles  only  asymptotically  approached 
those  of  fully  developed  turbulent  pipe  flow,  the  length  of  the  transi¬ 
tion  flow  regime  could  not  be  distinctly  defined. 

The  measured  dimensions  for  the  laminar  core  and  mixing  zone  are 
compared  with  theory  in  Figure  16.  The  dimensions  of  a  confined  jet  are 
quite  different  from  those  of  a  free  jet  (Reference  15).  In  the  present 
case,  the  dimension  of  the  laminar  core  is  a  function  of  the  initial 
velocity  ratio,  ;  consequently,  it  is  a  function  of  area  ratio 

There  exists  a  close  correlation  between  experiment  and  theory  of  these 
dimensions.  (Figure  16).  Also,  Figure  16  shows  close  agreement  between 
experiment  and  theory  for  the  boundary  of  the  mixing  zone . 

Measured  potential  velocity  is  compared  with  theory  in  Figure  11. 

The  experimental  results  for  the  growth  of  the  wall  boundary  layer  are 
compared  with  theory  in  Figure  15.  In  these  two  figures,  the  axial 
distance  was  converted  into  the  nondimensionalized  product  of  (  -i hr  ) 
ana  Reynolds  number  R,  where  K  was  based  on  the  radius  of  the 

ejector  Re,  and  the  induced  secondary  flow  velocity  .  In  Figure  11, 
the  potential  velocity  was  plotted  as  the  ratio  of  the  difference  between 
the  local  potential  velocity.  Up,  and  the  induced  secondary  velocity,  U^, 
to  the  induced  secondary  velocity,  U^. 

The  velocity  profiles  at  different  axial  locations  are  shown  in 
Figures  20-2&  A  comparison  of  the  theoretical  and  experimental  velocity 
profiles  forJA'  «  29.6  are  shown  in  Figure  17  for  the  modified  free- 
mixing  flow  regime.  Unlike  the  case  of  a  free  jet  flow,  the  value  or  n 
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approached  5.0,  as  can  be  seen  in  Figure  17e. 
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The  centerline  velocities  for  different  area  ratios  are  shown 
dimensionally  in  Figure  25.  They  are  shown  nondimensionally  and  compared 
with  the  centerline  velocities  for  free  jet  flow  in  Figure  26.  The 
centerline  velocity  for  free  jet  flow  is  inversely  proportional  to  the 
axial  distance.  (Reference  21).  The  inverse  of  the  centerline  velocity 
ratio,  ,  plotted  against  the  axial  distance  ratio,  -g-  ,  shows  a 
linear  variation.  (Figure  27).  For  the  confined  jets  in  this  report, 
the  experimental  value  of  these  quantities  could  be  approximated  by 
straight  lines.  The  formulae  representing  the  straight  lines  are  shown 
in  Figure  27.  The  slope  of  these  lines,  denoted  by  k,  increases  with 
area  ratio,  p.  ,  as  shown  in  Figure  28.  ks p.  approaches  infinity  k  tends 
toward  0.1362,  which  is  the  corresponding  slope  for  a  free  jet.  There¬ 
fore,  the  free  jet  flow  might  be  interpreted  as  a  confined  jet  flow  with 
/JL  equal  to  infinity. 

The  velocity  profiles  at  the  outlets  of  the  ejectors  were  measured 
and  plotted  nondimensionally  in  Figure  18.  The  data  are  compared  with 
theory.  It  can  be  seen  that  the  m  values  varied  proportionally  to  the 
length  ratio  L/D,  and  inversely  proportionally  to  the  area  ratio,  and 
that  the  m  value  varies  proportionally  to  the  ratio  of  the  length  ratio 
L/D  to  the  area  ratio  yU  ,  .  (Figure  19).  This  phenomenon  is 

explained  as  follows:  the  ratio,  ,  could  be  interpreted  as  an 

indication  of  relative  importance  of  the  frictional  forces  in  the  flow 
fi^ld;  i.e.,  skin  frictional  forces  per  unit  momentum  issued  from  the 
nczzle,  and  this  ratio  is  the  inverse  of  Reynolds  Number  based  on 
significant  variables  to  be  determined.  If  the  skin  friction  is  relative¬ 
ly  important,  the  velocity  profile  is  fuller  than  that  for  a  flow  field 
where  skin  friction  is  not  relatively  important.  The  uncertainty  raised 
herein  is  an  interesting  problem  for  future  research. 


Unlike  the  case  of  a  free  jet,  the  pressure  was  not  constant  in  the 
flow  field,  as  was  measured  and  plotted  in  Figure  29. 


B.  Conclusions : 

It  is  the  author's  belief  that  the  mixing  process  within  the  ejector 
is  an  interchanging  of  momentum  and  energy  among  each  stratum  of  the  flow 
field.  As  soon  as  the  ambient  air  is  entrained  into  the  ejector,  the 
interchanging  of  momentum  and  energy  occurs.  Particles  or  groups  of 
particles  in  the  flow  field  which  have  possessed  less  momentum  and  energy 
will  recieve  some  amount  of  external  momentum  and  energy,  and  on  the 
other  hand,  particles  or  groups  of  particles  which  are  at  higher  energy 
and  momentum  level  will  lose  some.  It  is  beyond  the  scope  of  this  report 
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Co  survey  exactly  the  causes  and  consequences  of  this  interchange  of 
energy  and  momentum. 

a.  The  Structure  of  the  Flow  Field  in  the  Inlet  Length  of  the 
Ejector. 

1.  The  jet  flow  emerges  at  the  exit  of  the  nozzle  as  a  constant 
velocity  profile,  and  the  secondary  Induced  flow  velocity  is  also 
constant  throughout  the  rest  of  the  cross  section  area  at  the  initial 
plane.  There  is  a  near  discontinuity  of  velocity  at  the  rim  of  the  exit 
of  the  nozzle.  This  near  discontinuity  in  the  velocity  profile  is 
smoothed  gradually  downstream. 

i.  A  laminar  core  starts  at  the  exit  of  the  nozzle,  and  its 
radius  decreases  linearly  with  respect  to  the  axial  coordinate  until  it 
becomes  a  point  on  the  longitudinal  axis  of  the  ejector.  The  length  of 
the  laminar  core  was  about  6  nozzle  diameters  for  the  present  experiment. 
Flow  within  the  core  is  essentially  uninfluenced  by  the  confinement  of 
the  Jet  flow.  The  laminar  core  in  the  region  of  the  turbulent-mixing 
flow  can  be  described  by  the  formula  derived  by  Helmbold. 

if  -  2.5(i-*0 

3.  The  Jet  flow  starts  mixing  with  the  surrounding  potential  flow 
at  the  exit  of  the  nozzle.  The  potential  flow  is  uninfluenced  by  the 
viscous  effects  of  the  walls  beyond  the  wall  boundary  layer.  The  mixing 
zone  spreads  outward  linearly  until  it  meets  the  boundary  layer  of  the 
ejector.  The  radius  of  the  mixing  zone  can  be  expressed  as 

*  =^i5(^olllz3  +  o.zs(^T-  0+  / . 

4.  The  wall  boundary  layer  of  the  ejector  is  infinitely  thin  at  the 
initial  plane  and  it  grows  until  meeting  the  boundary  of  the  turbulent 
mixing  zone.  The  growth  of  the  wall  boundary  layer  can  be  described  in 
terms  of  potential  flow  velocity 

5.  The  velocity  profile  within  the  wall  boundary  layer  of  the 

ejector  assumes  the  parabolic  shape  ~  , 

of>‘. 
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6.  Beyond  the  region  of  turbulent  mixing  the  potential  flow 
velocity  can  be  described  by  the  following  equation: 

=  ir  r  fCnJ 

7.  The  velocity  profile  within  the  turbulent-mixing  region  can  be 
expressed  by  Cornish's  extension  of  Cole's  two-dimensional  wake  law, 

^=i[i-cosn{>-(i-(asnj|}n'j  } 

where  n  varies  from  1.0  to  5.0  beyond  the  core  region  and  up  to  the 
location  where  the  boundary  of  the  turbulent  mixing  meets  the  edge  of 
the  boundary  layer  of  the  walls.  For  this  equation  to  be  applicable 
beyond  the  laminar  core  region,  w  has  to  be  changed  to  r  and  w^  to  r^. 


b.  Centerline  Velocity  of  the  Ejector. 

Centerline  velocity  is  inversely  proportional  to  the  axial 
distance  and  can  be  expressed  as 

u _ L-j _ 

Vo"  CONSTANT  +  -fc  (-§-) 

for  each  value  of  area  ratio. 


c.  Velocity  Profile  at  the  Outlet  of  the  Ejector. 


Velocity  profile  at  the  outlet  of  the  ejector  can  be  expressed 


as 


Ha-/*.) 

Uz  \Re)  , 

where  m  is  a  positive  real  number. 


d.  Velocity  Ratios  and  Augmentation  Factor. 

The  augmentation  factor,  Q  ,  and  the  velocity  ratios,  j 
andcX^  ,  car.  be  expressed  in  terms  of  area  ratio. 
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For  future  research,  it  would  be  interesting  to  explore  the 
characteristics  of  ar  ejector  with  different  primary  and  secondary 
fluids  with  different  temperature  ratios.  It  would  also  be  very  useful 
to  study  the  characteristics  of  this  arrangement  under  dynamic  conditions; 
i.e.,  with  the  secondary  fluid  flowing. 

For  simplification  of  design,  future  studies  of  the  jet  ejector 
flow  should  include  the  possibility  of  the  characteristics  of  the 
device  being  a  function  of  only  one  parameter. 
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Descriptive  Layout  of  Apparatus 
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Figure  3a.  Thrust  Force  Measuring  Unit. 


Figure  3b.  Jet  Flow  Unit,  Nozzle,  and  Supporting  Rig. 
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Figure  4.  Extension  for  the  Nozzle  Exit. 


Figure  5.  Pitot  Tube  Probe  and  Traversing  Mechanism. 
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Figure  6.  Velocity -Measuring  Setup. 
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Figure  7.  Thrust  Forces  of  Nozzles. 
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Figure  9.  Effect  of  the  Length  Ratio  on  Augmentation  Factor. 


Figure  10.  Effect  of  Gap  Ratio  on  Augmentation  Factor 
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Theoretical  Boundary  of  the  Turbulent  Mixing  Zone  and  Comparison  of 
Theoretical  and  Experimental  Thickness  of  Boundary  Layer  of  the  Walls 


Comparison  of  Theoretical  and  Experimental  Velocity  Profiles  in 


Figure  18. 


Comparison  of  Theoretical  and  Experimental  Velocity  Profiles  at  Outlet 
of  Ejector. 


Figure  19. 
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Figure  22.  Velocity  Profiles,  ytf  =  47.0. 
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Figure  26.  Comparison  of  Theoretical  and  Experimental  Centerline  Velocity  Ratio 
of  Ejectors  and  That  of  a  Free  Jet. 
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TABLE  1 

THE  DIMENSIONS  OF  THE  EXTENSIONS  OF  NOZZLES 


Nozzle  No. 

Inside  Diameter  of  Exits 

Extension  Lengths 

(inches) 

(inches) 

1 

0.304 

0 . 5*» ,  1.15,  2.12,  4.13 

2 

0.520 

0.60,  1.10,  2.10,  4.15 

TABLE 

THE  DIMENSIONS  OF 

2 

THE  EJECTORS 

Ejector  No. 

Inside  Diameter 

No.  of  Short 

Length 

(inches) 

Ejectors 

(inches) 

1 

1.04 

5 

11.40 

2 

2.08 

6 

15.60 

3 

3.12 

6 

17.15 

4 

2.82 

- 

42.00 

51 


TABLE  3 


AXIAL  LOCATIONS  WHERE  VELOCITY 
DISTRIBUTIONS  WERE  MEASURED 


Ejector  No.  Axial  Distance  from  Initial  Plane  in  Inches 


0.00, 

1.14,  2.14, 

3.14, 

5.14, 

7.14, 

9.24, 

11.14 

0.00, 

13.16 

1.16,  2.16, 

,  15.20 

3.16, 

5.16, 

7.16, 

9.16, 

11.16, 

0.00, 

1.00,  2.03, 

3.01, 

5.00, 

9.02, 

11.02 

,  13.00 

4  At  an  interval  of  0.25  inches  for  the  first  12 

■stations,  0.50  inches  for  the  next  28  stations, 
and  2.00  inches  for  the  last  6  stations. 


